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Abstract — We provide the solution for optimizing the power 
and resource allocation over block-fading relay-assisted broad- 
cast channels in order to maximize the long term average 
achievable rates region of the users. The problem formulation 
assumes regenerative (repetition coding) decode-and-forward 
(DF) relaying strategy, long-term average total transmitted power 
constraint, orthogonal multiplexing of the users messages within 
the channel blocks, possibility to use a direct transmission (DT) 
mode from the base station to the user terminal directly or 
a relaying (DF) transmission mode, and partial channel state 
information. We show that our optimization problem can be 
transformed into an equivalent "no-relaying" broadcast channel 
optimization problem with each actual user substituted by two 
virtual users having different channel qualities and multiplexing 
weights. The proposed power and resource allocation strategies 
are expressed in closed-form that can be applied practically 
in centralized relay-assisted wireless networks. Furthermore, 
we show by numerical examples that our scheme enlarges the 
achievable rates region significantly. 

Index Terms — wireless networks, relaying, resource allocation, 
multiuser diversity, long-term-evolution (LTE) 



I. Introduction 

Although the relay channel was investigated by the infor- 
mation theory researchers long time ago (JT], |]2]), the topic 
of cooperation/relaying schemes have recently become a very 
active research area within both the information theory as well 
as communications engineering societies. Few examples of 
recent works among many others are 0, H, A3, 0, [j?], 10, 
[|9l , and iflOl . It is well-understood now that relaying strategies 
can improve the coverage of wireless networks by providing 
higher data rates or better transmission reliability to terminals 
at the edge of a wireless cell (i.e. terminals which receive low 
signal power from the base station). 

Relaying technologies are also becoming part of the 
telecommunication standards ifTTI . Although many advanced 
schemes based on the cooperation of the mobile users to 
help each other are being studied in the literature, the first 
actual deployment step which will take place within the 3GPP3 
Long-Term Evolution (LTE)-Advanced (c.f. JE), O, El) 
standard is based on fixed access point^l to do the relaying 
and within a centralized scheme in which the e-nodeB (base 
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2 The relays will be base stations but without a wired or microwave 
connection to the backhaul network. 



station with backhaul connection) takes the scheduling and 
resource allocation decisions. One major objective in 3GPP 
evolution is to utilize the scarce wireless system resources ef- 
ficiently because achieving the high Quality-of-Service (QoS) 
targets through over-provisioning is uneconomical due to the 
relatively high cost for transmission capacity in cellular access 
networks [15|. 

Our objective here is to obtain the optimal (in information 
theoretic perspective) resource allocation schemes but with 
applying system constraints that are relevant to the LTE- 
Advanced standard so that it can be applied practically in 
the first introduction of relays to the wireless cellular systems 
industry. We have been able to derive the optimal power, 
resource allocation and scheduling polices that are provably 
based on closed-form formulations which are practical for 
implementation. Our proposed resource allocation schemes 
provide an integrated solution to exploit the cooperative (re- 
laying) diversity gains [3] as well as the multiuser diversity 
gains fl6l . 

Optimal dynamic resource allocation over fading channels 
has been investigated in the literature for non-cooperative 
(i.e. with no relaying) wireless systems in IfTTI for the single- 
user case, in |[T8l for multi-access (many-to-one) channels, and 
in lfl9ll and |20| for broadcast (one-to-many) channels. Optimal 
resource allocation for fading relay channels has been studied 
in lETl . Il22l for single source and single destination node 
case. Resource allocation for broadcast-relay channels was 
also treated in |23| and [24], where the problem setup assumes 
that users terminals cooperate together (i.e. act as relays to 
help other users). This is a different problem setup than the 
problem in this paper in which we assume that fixed-relays 
are used to assist in the transmission without being destination 
nodes themselves. In l25ll and 1261 . joint power and resource 
allocation over multiple access relay-assisted channels was 
studied for a constant channel realization scenario, where it 
was demonstrated that joint allocation of power and channel 
resources can enlarge the achievable rate region. 

In this paper, we consider the broadcast relay-assisted chan- 
nels under block-fading conditions (i.e. over many channel 
realizations). In E7] , ll28l . and ll29l relay-assisted broadcast 
(downlink) channels with centralized scheduling were consid- 
ered with different performance metrics such as throughput - 
guarantees or fairness measures. In our work, we tackle the 
problem from an information-theoretic perspective in which 
we aim to maximize the achievable rate region. This is equiv- 
alent to the problem of minimizing the transmitted power to 
achieve requested rate demands. To the best of our knowledge, 
maximizing the achievable rate region of block-fading relay- 
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assisted broadcast channels for the case of applying half- 
duplex regenerative decode-and-forward (RDF) J3] relaying 
strategy and orthogonal multiplexing of user messages within 
the channel blocks has not been treated in the literature. As 
well-known form the information theory literature, the achiev- 
able rate region can be enlarged (i.e. improved) by using non- 
orthogonal (superposition-based with successive interference 
cancellation at the receivers) transmission strategies as well 
as more advanced relaying strategies^ than RDF. However, 
we restrict our optimization problem with systems constraints 
which are favorable for a practical deployment. 

We formulate our optimization problem assuming par- 
tial channel-state-information (CSI). Similar to the definition 
adopted in |9|, by partial CSI we mean that the transmitter 
(source node) knows the channel state amplitudes only (i.e. the 
channel power gain), while the receiver knows both the chan- 
nel amplitude and phase. This can be usually accomplished 
by separate low rate feedback channels. According to LTE 
specifications [13], CSI is obtained by various means and 
forwarded to the base station, where the resource allocation 
decisions are done. The scheduling information are provided 
to the other nodes in the network using dedicated control 
channels. Furthermore, we apply a single long-term average 
total (i.e. sum) power constraint in the formulation of the opti- 
mization problem instead of using separate power constraints 
for every involved node (the source and relay nodes). This 
approach has been adopted in some works in the literature 
such as (9). The main reasons for applying a single power 
constraint are: (i) using separate power constraints for every 
node is less dynamic than using a sum-power constraint for all 
nodes since the solution of the latter formulation involves the 
optimal power distribution among the nodes to maximize the 
required objective (maximizing the achievable rate region in 
our case). On the other hand, with separate power constraints, 
the power distribution among the nodes is fixed beforehand 
and, as a result, we lose one factor that can help us to 
maximize our objective, (ii) we believe that the sum-power 
constraint is relevant in practice since the relay nodes are 
fixed access points that are not limited in energy supply (i.e. 
they are not running on batteries like mobile handsets). So, 
we can distribute the power among the relays based on the 
users distribution in the cell, and hence we can involve one 
relay more than the others if there are more users in its 
vicinity. Note that if no relays are involved at all, the base 
station will transmit all the power anyway. In our problem 
formulation, we are just re-distributing some of that total 
power among the relays in the most efficient way, (iii) solving 
the optimization problem with many power constraints will 
be based on bisection-based methods over many dimensions 
similar to the solution in ll25l for the multiple-access case, 
while using a sum power constraint will result in closed-form 

3 For example, full-duplex relaying strategies, which require the relay 
terminal to be able to receive and transmit simultaneously on the same chan- 
nel, gives better throughput gains than half-duplex strategies. However, this 
requirement is difficult to implement in practice because of the complexity of 
providing electrical isolation between the transmitter and receiver circuitries. 
Thus, half duplex relaying strategies are favorable from practical point-of- 
view although they result in multiplexing loss because the same message is 
transmitted over two orthogonal slots. 



solutions which are simpler and more practical to be applied 
in real systems. 

Following this introduction section, we present in Section HIl 
the channel model and the optimization problem formulation 
as well as the applied mathematical notation. In Section [ill] 
we provide the solution of the optimization problem. Despite 
the relative large number of optimization variables involved in 
the problem, we demonstrate that the problem can be solved 
by: (i) characterizing the maximum achievable rates using 
RDF links with optimal power allocation over the source 
and the relay, and (ii) transforming the problem into an 
equivalent "no-relaying" broadcast channel with each "actual" 
user replaced by two "virtual" users having different channel 
qualities and multiplexing weights. In order to understand 
the mathematical solution steps, the reader is advised to go 
through [19 Section III-B] in which the orthogonal "non- 
cooperative" broadcast channel was considered. We give the 
solution for our problem which includes the closed-form 
policies to select the best relay, schedule the users across 
the resource units, choose the optimal transmission mode and 
control the transmission power. We include also the case of 
optimizing resource allocation over each channel realization 
separately and use it for comparison with the optimal case. 
We provide numerical examples in Section [IV] to demonstrate 
the advantages of our suggested resource allocation scheme. 
Then, we summarize the main conclusions of our work in 
Section [Vj 

II. System Model and Problem Formulation 
A. Channel Model 



Orthogonal Regenerative 1 
Decode-and-Forward Link 



Direct Transmission Link 



Fig. 1. System diagram of relay-assisted broadcast channels with one 
source node, L relay nodes and M destination nodes. The message sent to a 
destination node can go through half-duplex RDF link with the assistance of 
one of the relay nodes, or through a direct transmission link from the source 
node to the destination node without the assistance of the relays. 



We consider the Af-user relay-assisted broadcast system 
shown in Fig. Q] The system consists of one source node 
(S), M destination nodes (D), and L relay nodes (i?) which 
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can be used to assist in the transmission from the S-node 
to any of the D-nodes. The i?-nodes are connected to other 
nodes by wireless links only. Any message transmitted by 
the S-node is destined to only one £>-node. We assume that 
the message transfer from the S'-node to a D-node can go 
through two possible modes: (i) a direct transmission (DT) link 
without the assistance of any of the i?-nodes, (ii) a half-duplex 
regenerative (i.e. repetition coding) decode-and-forward (DF) 
link in which one of the relays assists in the communication 
between the source and the destination. The S'-node represents 
the base station, while the D-nodes correspond to the mobile 
users' terminals. 

We assume a block-fading channel model in which the 
air-link resource grid is divided in both time and frequency 
domains into small blocks called resource units (RUsjj. The 
channel is assumed to be constant within one RU, but varying 
(fading) independently across the RUs in the air-link grid. The 
resource units are orthogonal to each other (non-overlapping). 
Furthermore, we assume that all messages are multiplexed 
orthogonally in all RUs. Irrespective of their transmission 
mode (i.e. DF or DT), one or many users can receive in the 
same RU (subject to optimization). However, if more than 
one user is receiving messages in the same resource unit, 
the messages of the users are multiplexed orthogonally by 
frequency division^ A DF link is divided into two sub-units 
occupying the same frequency band but having orthogonal 
time division multiplexing. Fig. [2] shows how the air-link 
resource grid is divided, and an example of how the users 
may be scheduled across the RUs. One of the possibilities is 
that two independent messages are sent to the same D-node in 
the same RU, but with different transmission modes; i.e. one 
message is sent through a DT link and the other one using a 
DF link. Several other possibilities are illustrated in Fig. [2] 

Additionally, we assume that the channel realizations 
(blocks) are large enough (i.e. slow fading assumption) so that 
each codeword destined to any user can be transmitted over 
one channel block with close-to-capacity limits rate. 

B. Problem Formulation 

In this work, we provide answers to several fundamental 
questional concerning the resource allocation schemes to be 
applied in the system: 

1) When (i.e. under what conditions) is a DT link optimal? 
When is a DF link optimal? When is the orthogonal 
multiplexing of both modes optimal? 

4 As well-known from LTE specifications, the time-frequency resource grid 
is divided into resource units (RU) which are allocated flexibly to user 
terminals [14|. A pair of adjacent resource units (0.5 ms each) is allocated to 
one terminal based on the scheduler decisions. 

5 As known from LTE specifications, one RU spans several OFDM subcarri- 
ers. So, in principle, we can let multiple users share the same RU orthogonally. 
However, according to LTE specifications, one RU can be used by one user 
only. So, it may appear that our problem formulation is not practical to be 
applied in LTE systems. However, as discussed in Section IfTfl the solution of 
the optimization problem with the orthogonality constraint involves that only 
one D-node receives in one channel block (RU unit) using either DT or DF. 

s Some of these questions may have been discussed in the literature for 
other channel models or problem setup. We are interested here in the optimal 
resource allocation scheme within the assumptions and channel model adopted 
in this paper. 



Time 



R 1 — > D 1 
S — > R 1 &D 1 


S — > D 


S — > D 3 


R i ~ _> D 3 


R 2 " -> D 2 


R 2 _> °1 


S -- > D 


S— >R 1 &D 3 


S — > R 2 &Dj 


S — > R 2 &D ) 


Resource Unit 
"Channel Block" 


R 1 ~ _> °S 


S — > R^Dg 



Fig. 2. Dividing the total air-link resource grid into small blocks called 
resource units, which are assumed to be flat faded (i.e. the bandwidth of a 
RU and its time duration are smaller than the coherence bandwidth and the 
coherence time of the fading channels). Within one RU, one or more users can 
receive a message from the source node. All transmissions are orthogonal to 
each other. Based on the channel conditions, the messages and their associated 
transmission modes can be scheduled flexibly to the users terminals. 



2) How should we select the relay that is most capable of 
supporting the transmission to a D-node? What is the 
optimal policy to allocate the power for the S-node and 
the i?-node over a DF link? 

3) How should the scheduler process the CSI to determine 
the number of messages to be transmitted in a RU and 
their transmission modes? What is the optimal policy 
to allocate the power and the channel ratios for every 
message? 

Furthermore, in a flexible system where the users can 
have different requested rates dependent on, for example, the 
supported services, the resource allocation schemes should be 
flexibly adjustable to operate at any of the possible operating 
points of the system. In this paper, we solve the resource 
allocation problem to operate at any pre-determined point 
on the achievable rates region. The achievable rates region 
is defined as the set of all long-term average rate vectors^] 
R = [_Ri R.2 ■ ■ ■ Rm] achieved by the nodes D\,D2, ■ ■ ■ , Dm 
such that the long-term average sum (of all nodes) power 
density constraint P is not exceeded. To simplify the math- 
ematical formulation of the problem, we assume that B and 
P represent spectral densities (i.e. bits/sec/Hz and Joul/sec/Hz 
respectively). 

The optimum points within the achievable rates region are 
those that are located on the boundary surface. The latter can 
be characterized as the closure of the parametrically defined 

7 We use boldface to indicate vectors. 
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surface 



{R(/x):/xe^,5> 4 = l} 



(1) 



where for every weighting-factors vector /x, the rate vector 
R(/x) can be obtained by solving the optimization problem: 



M K 

max'y^iiiRi, subject to — P[k] < P 



(2) 



fc=i 



where K — > oo is the total number of channel blocks, 
and M is the number of active users. Ri is the long-term 
average (i.e. averaged over all channel blocks) rate of user 
i. The relationship between R and P will be obtained in 
Section [HI] For simplicity, we use a single index fc to refer 
to a RU although the resource grid is divided in both time 
and frequency domains. We assume without-loss-of-generality 
that all channel blocks (RUs) have identical frequency band- 
width and time duration. Pi[k] and P[k] are respectively the 
sum (of the source node and the relay node) power density 
(Joul/sec/Hz) used specifically to transmit to node Di, and 
the sum power density (Joul/sec/Hz) transmitted (including all 
receiving nodes) during channel block fc. 

All power and resource allocation polices proposed in this 
paper are presented as functions of the weighting factors vector 
because it defines the specific operating point of the system. 
The selection of /x to meet the constraints of the provided 
services (applications) is a different topic that is not discussed 
in this worlo Few examples of the many possible approaches 
suggested in the literature to select the specific operating 
point of the system are (i) the fairness-based approach, such 
as the proportional fairness scheduler 11321 and the flexible 
resource-sharing constraints scheduler [1331 . (ii) the utility- 
maximization-based approach 11341 . and (iii) the Quality-of- 
Service (QoS) constraints based approach [35|[36|. 

C. Mathematical Notation 

The notation Ri [fc] means the achievable rate (bits/sec/Hz) 
at node Di during one RU, which has index fc. In our mathe- 
matical notations, we use superscripts to differentiate between 
DT and DF links. Similarly, we use subscripts to indicate the 
nodes involved. We denote the relay node associated with 



destination node D, as R 



Furthermore, we denote the 



channel access ratios (i.e. bandwidth ratio) of user Di within 
channel block fc as Tf T [fc] for the DT link, and rf F [fc] for the 
DF link. It should be clear that T,iLi(T t pT [k]+T? F [k]) = 1. 
The channel gain between two nodes is represented by h. For 
example, h s di is the channel quality between the source node 
and the destination node (Di). Fig. [3] shows an example of 
one RU (with index fc) with a DF link to one node Di as well 
as a DT link to another node Dj multiplexed orthogonally in 
the frequency domain. The relations between the transmitted 
signals (x DF [k), x® T [k\) and the signals (y[k\) received by 

8 As discussed in 1301 and 1311 . there is no contradiction between the two 
objectives of (i) efficient resource allocation by designing resource allocation 
schemes leading to operating at the points on the boundary of the achievable 
rates region, and (ii) achieving fairness among the users as well as maintaining 
the QoS requirements, which can be done by controlling the operating point 
of the system based on proper selection of fj.. 



the destination and relay node associated with node Di in 
each of the orthogonal sub-channels within the resource unit 
fc are shown in Fig. [5] Here, zf F [fc] and z^ F ^ [fc] represent 
the additive zero-mean white circular complex Gaussian noise 
with variance a 2 at Dj and Rm) respectively. 
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Fig. 3. An example of the orthogonal division of one resource unit k to 
transmit to users Di (using DF mode with assistance of node i?.;(ij) and Dj 
(using DT mode). The received signals at the associated nodes are shown as 
functions of the transmitted signals. 



We use the notation 7 = U- for the effective power gain of 
a given channel, where N a is the noise power spectral density. 
We assume that the fading processes of the channel gains 
("fsdi, 7sr !(i )' lr l(i) di) are independent of each other, stationary 
and have continuous probability density functions, / 7 (x). In 
the numerical examples throughout the paper, we assume the 
fading processes have Rayleigh@ or Rice0 distributions l37l . 



III. Solution Structure 

Looking back at our optimization problem (f2), the achiev- 
able rate Ri [fc] is the sum of the rates achieved by the DT and 
DF links multiplied by their relative channel access ratios: 

n[k]Ri[k] = rP F [k]Rf F [k] + r, DT [fc]i?f T [fc] (3) 

The maximum possible achievable rate (bits/sec/Hz) of user 
i through DT link is given by: 



R DT [k] = log (1 +7.dJfc]^H fc ]) 



(4) 



for additive white Gaussian receiver noise, where (0) is the 
Shannon capacity for the AWGN channel. With adaptive mod- 
ulation and coding, a rate close to capacity can be achieved 
(e.g. [ 38 1 ). In practice, wireless systems support a set of 

9 /^ ayleigh (x) = i exp a; > 0, 7 is the average effective power 

gain of the channel. 

10 /^ c (z) = ^exp (-« - iS±l*) I f 2y /^^\ x > 0, 7 is 

the average effective power gain of the channel, and k is the ratio of the 
power received through the line-of-sight path to the power received through 
the non-light-of-sight path. 
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discrete rate values rather than a continuous range. However, 
we use the idealization to relate "power" and "rate"; the 
relative performances^ will carry over into practice. 

The achievable rate by regenerative decode-and-forward 
relaying is known from the literature (e.g. [f3]): 

R? F [k] = 

j ±log(l + 7sn(i Jfc]P^f[fc]), 
mm \ I log (l + ladi [k]P s D d F [k] + lrmdi [k]P% F di [k]) 

(5) 

The first term in (0 is the achievable rate by the relay 
node in the first time slot. Since, in RDF, the relay has to 
decode the signal in order to retransmit it in the second time 
slot, the achievable rate by the destination node is upper 
bounded by the "decodable" rate limit at the relay. The 1/2 
factor which is multiplying the log function is due to the fact 
that the information is conveyed to the relay in half the total 
time assigned for transmitting the signal to the _D-node. The 
second term in (O is the maximum achievable rate by the D- 
node using Maximal-Ratio-Combining (MRC), e.g. fl6| , of 
the signal received from the source in the first slot, and the 
signal received from the relay in the second time slot. 

In the sequel, we assume a single relay selection per-user 
and per-resource unit in RDF transmission mode. This is 
appropriate from a practical point of view especially that we 
assume partial CSI (no phase information) available at the 
transmitters. In Appendix[A]we extend the solution to the case 
of multiple-relay selection per-user which requires full CSI 
(i.e. including phase information) in order to enable coherent 
(in phase) transmission of multiple relays. 

The set of optimization variables for the problem (0 
includes the power allocated to each node as well as the 
corresponding channel access ratios for every transmission 
mode in every channel block k. We need first to characterize 
the achievable rate by RDF links in order to solve (|2). 

A. Achievable Rate by Regenerative Decode and Forward 

Our objective here is to maximize R F>F [k] §5§ with respect 
to the sum power allocated to the DF link: 



.) ■ " •)-■</, " i" j (6) 

The 1/2 factors in © are used because we present P as 
(Joul/sec/Hz), and the source and relay nodes are transmitting 
in half the total time. For simplicity, we will drop the user in- 
dex i as well as the channel block index fc in our mathematical 
formulations in this sub-section. 

We first define when the DF link can be useful. By the 
notion "useful DF link" we mean that it could (for some power 
range) support higher data rate than a DT link given that both 
transmission modes are allocated the same total power. 

Proposition 1: If 7 sr or j rd is less than ~f sd , then transmit- 
ting using a DT link supports higher data rates than using a 

1 1 The "absolute" performance of a combination of specific modulation and 
coding schemes can often be approximated by l|4j as well. An "acceptable" 
residual bit or frame error rate will often be achieved by a practical scheme 
with some (fairly constant) power-offset against the theoretical "zero-error" 
curve given by l|4j- 



DF link, which means that the DF link is NOT useful in this 
case. 

Proof: If j sr < j sd , then the relay node can decode at 
rates less than the rates achievable by the DT link. Similarly, 
if Ird < 7sd, then allocating the power in the second slot 
to the source node to retransmit the codeword is better than 
allocating the power to the relay node. ■ 
The next step is to solve the problem: 

max R DF -mm( ? bg i 1 + 7s ^ ^ ' 

(7a) 



subject to 



Isr > Isd, Ird > Isd, 



DF 



DF 



= P 



DF 



(7b) 



Proposition 2: The optimal allocation of the source power 
and the relay power over a useful DF link can be obtained by 
making the two terms in dTat equal each other. 

Proof: The first term in (fTa) (i.e. \ log (l + j sr P^ F )) 
is a monotonically increasing function of P S D . 
On the other hand, the second term in (|7a] i 
(i.e. \ log (l + 7 sd P® F + ^ rd P^ F )) is a monotonically 
decreasing function of P® F because 7^ > -f sd , and the 
sum of P® F and P F>F is constant (equals 2P DF ). Thus, to 
maximize the minimum of the two terms in (l7ai i. we should 
make them equal. ■ 

Hence, the power allocation over a useful DF link should 
be as follows: 



P. 



DF 



2P 



DF 



DF 



2P 



DF 



P. 



DF 



Jrd 



where j' sr and j rd in (© are defined as: 



7s 



Isr , . Ird 
, Ird = 

Isd Isd 



(8) 



(9) 



Furthermore, the achievable rate over a useful DF link 

(i.e. Ysr > 1 and ^' rd > 1) is: 



1 



where 



R ut {P ut ) = ^og(l + 2a lsd P 



Isrlrd 



DF\ 



(10) 



(11) 



l'sr + Ird - 1 

The parameter a is actually the power gain that the RDF 
link is capable of providing. It should be clear from (ITTt that 
if a relay is useful, then a > 1. The result in ( TTOb is interesting 
because it fits with the expectation that a half-duplex relaying 
strategy provides a power gain (i.e. beamforming gain) at 
the cost of a loss (by half) in the degrees-of-freedorr[3 
(i.e. multiplexing gain). The beamforming gain is obtained 
because the receiver decodes the signal at higher Signal-to- 
noise-ratio (SNR) although the same total power is used for 
the transmission. The loss in the multiplexing gain is a result 
of transmitting the message over two time slots. Thus, RDF, 
and all half-duplex relaying strategies in general, can be useful 

12 Readers who are not familiar with the notions of degrees-of-freedom, 
bandwidth-limited and power-limited regions of the AWGN channel capacity 
formula can refer to [16. Chapter 5]. 
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for terminals which are operating in the power limited region 
of the channel capacity where the achievable rate (bits/sec) 
has almost a linear relation with SNR. This means that RDF 
can support higher rates than a DT link only for users with 
low SNR (i.e. close or at the cell edge). On the other hand, 
terminals operating in the bandwidth-limited region of the 
channel capacity (i.e. at high SNR) will be affected by the 
loss in degrees-of-freedom of RDF, and hence a DT link can 
support higher data rate in this case. Actually, it can be shown 
that other half-duplex relaying strategies such as amplify-and- 
forward (AF) have achievable rate performance similar to (TT~0T >. 
but with more complicated formulas to express the power gain 
a in terms of the channel qualities between the source, relay 
and destination nodes. RDF has the nice property that a is 
a function of the channel qualities only and independent of 
the total power allocated to the link. This is not the case with 
other relaying strategies such as AF. 

Now, we can define the criteria to select the unique relay 
that should be selected in the RDF link. It is the one that 
provides the best power gain in order to maximize the total 
achieved rate at the destination node given the total power 
allocated to the DF link. The best relay to be associated with 
node Di is Ri(i) where is defined as: 



= argmaxo! r ,(j. 

3 



7/ rj > 1 and ^ r ] di > 1 (12) 



In the remaining part of this paper, we will use (fTOt to 
characterize the performance of RDF links. We assume that 
proper relay selection is used (Q~2j, and based on the allocated 
power to the RDF link (P DF ), the source and relay power 
should be allocated according to (0. For simplicity, we will 
drop the index of the relay and just use ati[k] to indicate 
the best power gain for the RDF link. If none of the relays 
is useful, we will use the DT link only for transmission. 
Otherwise, an orthogonal multiplexing between the DT link 
and the DF link should be used (subject to optimization). 



where the virtual user j that is related to one of the transmis- 
sion modes of the actual user i is characterized by: 



Mi 

ILl 

2 



if DT 
if DF 



Vj 



(14) 



7,* if DT 
2j sdt ai ifDF 

where fii is defined in the problem formulation (f2), and 
a can be obtained using (TlTT i and (TT2l) . As well known 
|[T9l , if the cumulative-density-function (CDF) of the fading 
process is a continuous function (this is true in reality such 
as in Rayleigh and Rice fading conditions), then the optimal 
resource allocation policy is unique and at most a single 
user only is scheduled in each fading state (i.e. channel 
block). Furthermore, the power allocated to the scheduled user 
follows a water-filling approach. 

Thus, every RU should be allocated to only one user with 
only one transmission mode (i.e. either DF or DT and not 
both). The index m of the only virtual user that should be 
scheduled in channel block k (i.e. Ti[k] = 1 if I = m and 
n [k] = if I m, I = 1, 2, • • • , 2M) is: 



m = argmax(/ J -(P j [A;]) - X G Pj[k}) 



(15) 



where fj is defined in (l3[ . and Xg is the "power price' 
which should be controlled to maintain the average power 
constraint P. Pj [k] in ( TT3T > is dependent on the transmission 
mode of the corresponding actual user i related to the virtual 
user j: 

w, 1 



i + 



■j 



vj[k] 



(16) 



where (x + = max(a;, 0)). Once we obtain m according to 
(TOl l. we can determine the actual user that should be scheduled 
and the optimal transmission mode. 

Expressions (Tl2l . (TOT l. ( fT~6b and ^ provide valuable closed- 
form policies to select the best relay, schedule the users across 
the RUs, choose the optimal transmission mode and control 
the transmission power. 



B. Optimal Resource Allocation - Broadcast Channel with 
2M Virtual Users 

The results in Section IIII-AI leads to a very interesting con- 
sequence that we can transform our problem into a broadcast 
(with no relays) channel with each actual user replaced by two 
virtual users having different channel qualities and multiplex- 
ing weights, corresponding to the two different transmission 
modes (i.e. DT or DF) of the actual user. The only difference 
between our problem and the original broadcast channel is 
that we will have two classes of users in terms of the relation 
between the achievable rate and the allocated power. However, 
this does not change the structure of the problem solution and 
the interesting closed-form solution of it lfl9l . The achievable 
rates multiplied by the corresponding multiplexing factor of 
the 2M virtual usersf^l will have the form 

f j (P j )=cj j log(l+r ]j P j ) (13) 

13 If some users do not have useful relay links, then the total number of 
virtual users will be less than 2M. 



C. Optimal Resource Allocation with Constant Power per 
Channel Block 

If the resource allocation is optimized for each channel 
block independently from other channel blocks and irrespec- 
tive of the scheduled users or selected transmission mode in 
each channel block, the corresponding optimization problem 
is: 

2M 

max y^rAklfjiPAk]), (17a) 

{r J [k],P,[k]j=l,...,2ilf}£j 31 Un Jl lh 

subject to 

2 M 

^2Tj[k]Pj[k] < P[k]. (17b) 

i=i 

where M is the number of active users, and fj is the achievable 
rate multiplied by the weighting index as defined in dT3l . Note 

14 The optimal value of Xq in order to achieve a desired average power 
level is dependent on the channel statistics. One option is to control in 
real-time based on actual channel measurements |35|. 



7 



that since we assumed in our problem definition in Section lTl-BI 
that R and P represent spectral densities (i.e. their units are 
bits/sec/Hz and Joul/sec/Hz respectively), the actual short-term 
average rate and power (averaged within one channel block 
k) should be obtained by multiplying with the corresponding 
channel access ratio (within channel block k) Tj[k] for each 
virtual user j. 

The solution of (T% can be obtained by applying the concept 
of virtual users over a no-relaying broadcast channel that is 
discussed in Section IIII-BI The solution of the problem lfl9l 
is that either one or two virtual users (could be related to the 
same actual user or to two different actual users) are scheduled 
within the channel block k. 

Depending on the value of P[k], the solution of the general 
case of problem (TTTb is either one or two elements in P 
and r are non-zero. The domain [0, oo) of the functions fj, 
j = 1, • • • , 2M will be divided into adjacent intervals which 
are alternating between (i) intervals in which the solution of 
the optimization problem is that the optimization variables 
of only one function (i.e. virtual user) should be non-zero, 
given that P[k] belongs to the corresponding interval, and (ii) 
intervals in which the optimization variables of two functions 
are non-zero, with the values of the P variables of these 
two functions identical to the values of the two end-points 
of the corresponding interval, and the r variables for these 
two functions dependent on the relative location of P[k] with 
respect to the interval end-points. An algorithm was provided 
in |[T9l pp. 1089, first column] to obtain the solution. 

Fig. |4] provides an illustration of the solution to problem 
([17. . The example presented in Fig. |4]is for three virtual users 
(two actual users with one of them having no useful DF link). 
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Fig. 4. An illustrative example of problem 1171 with three virtual users. The 
dependence of the solution on the specific interval, to which P belongs, is 
illustrated. 

In addition to the optimal solution, we suggest also a sub- 
optimal solution which is simple and close to optimality. The 
near-optimal solution is: 

T m = 1, P m = P[k], where m — argmax/i(P[fc]) (18) 



This near-optimal solution is identical to the optimal solu- 
tion when only one virtual user is optimizing the problem. 
However, if two virtual users are involved in the optimal 
solution, the sub-optimal solution still gives near-optimal 
performance because the difference between the tangent line 
between the two functions, which optimize the problem, and 
the maximum of the two functions is usually very small. 

IV. Numerical Examples 




10" 1 10° 10 1 

Average SNR at receiver through DT Link 



Fig. 5. Comparisons of the achievable rates using (i) direct transmission only, 
(ii) orthogonal multiplexing of RDF and DT with optimal power allocation. 
The solid lines correspond to the case of constant total power per channel 
realization, and the dashed lines correspond to the case of optimal power 
allocation over all channel realizations. 

We provide in Fig. [5] numerical comparisons of the maxi- 
mum achievable rates in a single-user case with and without 
applying the optimal power allocation strategy studied in this 
paper. The simulation results were obtained for the scenario 
when the power gain over the source-destination link is 
Rayleigh faded, while the power gain over the source-relay and 
relay-destination is Rician faded with better average channel 
qualities than over the direct link. Such scenario is relevant 
when there is no line-of-sight (LOS) path between the source 
and the destination nodes, while the relay node is in a position 
where it has LOS paths to both of them. The simulation results 
in Fig. [5] were done for two cases: (i) 2^ = 5, ^ = 3, and 
(ii) 3^ = 10, ^ = 5. In the simulations, the used values for 

' ~tsd ltd 

the ratios between the power of the LOS path over the non 
LOS paths in the Rician faded channels are: n sr = 10 and 

Fig. [5] shows that potential gains can be obtained with the 
assistance of the relay. The gain in the achievable rate is 
higher at lower SNR, but still useful at mid to high SNR. 
Furthermore, we observe that using optimal power allocation 
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Fig. 6. Comparisons of the maximum achievable rates with relay assistance 
for two cases: (i) optimal power allocation, (ii) equal power allocation for the 
source and the relay node. 

Another important observation is that with equal power 
allocation between the source and the relay nodes, which is 
shown in Fig. [6] there is a significant degradation from the 
achievable rates with the assistance of the relay. Furthermore, 
it is worse than just using direct transmission when the SNR 
is mid to high range. This demonstrates the advantage of 
the proposed power allocation strategy because it switches 
between DF and DT based on the channel conditions, and 
thus it is always better than just DT over all SNR regions. 

The results in Fig. [5] and Fig. [6] were obtained assuming 
large number of channel blocks (realizations) in order to obtain 
good estimate of the average performance of the system. Fig. [7] 
shows the ratio of RUs - with respect to the total number of 
RUs used in the simulation - which have a relay link, a direct 
transmission link, or none for the case of (3^ = 5, ^ = 3) in 
Fig. |5J The results in Fig. |7J include the global power control 
case as well as the constant power per channel block case. 
It is demonstrated that the DF link is more used at low SNR 
and decreases gradually as SNR improves, while at high SNR, 
the DT link is more capable of providing higher throughput 
to the D-node. The "no transmission" case appears when the 
channel condition is so bad that the optimal power control 
strategy is not to transmit during such deep fading conditions. 
It is obvious that such case appears more frequently at low 
SNR. 

Fig. |8] displays numerical comparisons of the achievable 
rates region in a two-user case with and without the assistance 
of a relay for the two cases of optimal power allocation 
over all channel blocks, and constant power per channel 



1.5 




0.5 1 15 2 25 3 

First User Average Rate (bits/sec/Hz) 



Fig. 8. Comparisons of the achievable rates in a two-user case with and 
without the assistance of a relay, and with and without power control over 
the channel blocks. The links parameters are: -y S r = 10 (Rician faded with 
K sr = 10), *f 3 d 1 = 10 (Rayleigh faded), 7 r d 1 = 2 (Rician faded with 
K r d 1 = 2), 7 S( 2 2 = 1 (Rayleigh faded), 7 r d 2 = ^ (Rician faded with 
K rd2 = 5)- 



block. The simulation results were obtained using the channel 
quality parameters shown in the Figure caption. The results 
in Fig. [8] demonstrates that the RDF strategy improves the 
performance of the system, especially for users who have bad 
connection with the base station. Furthermore, we observe that 
using optimal power allocation over all channel realizations is 
advantageous, especially at low SNR. However, using constant 
power allocation per channel block has good performance at 
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higher SNR. 



V. Conclusion 



We have addressed the problem of optimal resource al- 
location over block-fading relay-assisted broadcast channels 
under practical system constraints including using orthogo- 
nal multiplexing of transmitted messages, and regenerative 
decode-and-forward relaying. We have formulated the opti- 
mization problem using a sum power constraint, which is 
a valid assumption in fixed-relays scenarios with no power 
limitations, and we have shown that this is sensible in order 
to achieve the maximum possible performance as well as 
to transform the problem into an equivalent "no-relaying" 
broadcast channel optimization which has simple closed-form 
solution that is practical to be applied in real-time systems. 
The optimal resource allocation strategy is to schedule only 
one user within one channel block (i.e. resource unit) and 
to transmit using either a direct transmission mode without 
the assistance of the relays or through a relaying link. The 
relay selection is based on the maximum power gain that 
can be achieved via relaying. The power allocated for each 
transmission follows a water-filling approach, and in case of 
a relaying link, a formula is derived to determine the power 
allocated to the source and the relay respectively. If constant 
power is allocated per channel block, the optimal resource 
allocation scheme is to schedule not more than two users 
within one channel block. Simulation results demonstrate that 
our proposed resource allocation scheme provides considerable 
throughput gains especially for users receiving low power from 
the base station. 

Appendix A 
Multiple Relay Selection with Coherent 
Transmission 

We extend the results presented in Section [III] into the case 
when multiple relays can be involved in the RDF transmission 
mode within the same resource unit (RU). In this case, all 
involved relays listen to the signal transmitted by the source 
in the first time slot, and then decode the signal and retransmit 
it coherently to the destination node in the second time slot 
such that the transmitted signals by the relays add "in-phase" 
at the receiver in order to maximize the received power at the 
destination Z3-node in the second time slot. This necessitates 
that in addition to the channel "power gain" information, the 
channels' "phase" information should be known as well at the 
transmitters. 

The multiple relay selection case can be viewed as if the 
relay nodes represent antennas of one "virtual" relay that has 
multiple distributed antennas. Thus, the relay-destination link 
in this case is a multiple input single output (MISO) channel. 
Furthermore, we can allow the source node as well to transmit 
in the second slot as one of the "virtual" antennas of the 
"virtual" multiple-antenna relay. The capacity and optimal 
resource allocation over MISO channel is well known in the 
literature (e.g. [ 16 1). The maximum achievable effective power 



i in channel block k is the sum of the effective power gains 
of the individual links from every antenna to the _D-node: 



7. 



DF MISO\u „, r 



jen*[fe] 



(19) 



where fij[fc] is the set of indices of involved relays in the 
transmission to Dj in channel block k. 

The maximum achievable effective power gain ( [T9l > can be 
achieved by in-phase transmission of the distributed antennas 
and by adjusting the power transmitted over every antenna j to 
be directly proportional to its associated effective power gain 
7 r d- Thus, the power (spectral density) allocated to the relays 
follows: 

Irjdi [k] 



P 



DF_MISO 



r-: di 



[k} = 



DF_MISO \U 



P, 



(20) 



where P^~^d* S ° [k] is the sum (of all antennas) power 

(spectral density) allocated to the second time slot of the 

RDF link, and l?f -^" S °W\ is obtained in (19}. Similarly, 

the power allocated to the source node in the second time slot 

of RDF link is pDF_MlSO\u\ _ ,„W*1 pDF MISO r 
oi Kur unit is r gd [K\ — DF _ M iso, k] -r r , . d 



O i [fc]° 



All relays of the MISO RDF link should be able to decode the 
transmitted signal by the source in the first time slot. Thus, 
the achievable rate region of the RDF MISO link is: 



DF MISO 



R 



mm < 




I _|_ jDF_MISO rn p/)/'_.W/N<;[ 



1 



1 FA-1 P DF 

/sr*n. rfci win L J sdi 

Mp DF_MISO [k y 



(21) 



%l DF_MISO J ;„j pDF_MISO 



Isdi 

^r n . [k] di L A, J- t r sii [k]di 

where 7,E F - A f [k] is the minimum effective power gain 
of the source-relay channels (among the involved relays in 
the MISO relay link). Since all relays should decode the 
information in the first time slot, the achievable rate in the 
first time slot is bounded by the worst channel condition 
among all relays. Our objective is to find the optimal power 
allocation in the first and second time slots of the RDF links, 
and additionally to select the relays to be involved in the 
distributed MISO relay link. We first observe the following 
facts: 

• All relays that have 7 sr < -f s d should not be involved 
in the MISO RDF link in order for the RDF link to 
support higher achievable rates than in the case of direct 
transmission (DT). This fact can be easily proven similar 
to Proposition 1 in Section IIII-AI 

• Proposition 2 in Section IIII-AI is also valid in the case 
of multiple relays selection. This means that, for a given 
set (fij) of selected relays for the MISO RDF link, the 
power allocated to the first time slot (P^ F - MISO ) and 
for the second time slot (P^j 1 ^ 1130 ) respectively follow 
(© and (O with the replacement of ^ sr and j r d by 
7^-^in and l?n di IISO respectively. Similarly, the 
achievable rate can be obtained using (fT0T > and (fTTT i. 



gain 7, 



DF_MISO 

r siilk]di 



[k] of the MISO channel to destination node The power gain a that the MISO RDF link can support can 
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be obtained, similar to ( fTTT i. by: 



, DF_MISO 



,,DF_M 1 SO 



(22) 



- 1 



We aim to obtain the set of relays that maximize d22t . 
It is obvious that a increases (or decreases) by increasing 
(or decreasing) 7^^- A ^f° or 7,^ F df //5 °- We can increase 
7irn _7 min by removing the relay with worst channel quality 
with the source. However, by removing this relay from the set 
of relays, we decrease lrsfdf IIS ° which is a sum of all relay- 
destination channels gains. Thus, we can adopt the following 
procedure to find the set of relays that produce the highest 
possible power gain ( 1221 : 

1) Start with the set of all relays that satisfy the condition 
Isr > 7sd> an d compute the achievable power gain d22t . 

2) Modify the set of relays by removing the relay that has 
the least j sr , and re-compute the power gain. 

3) Repeat step 2 until you get a set of one relay only. 

4) Compare the power gains of the examined sets of relays' 
selection and choose the one that support the highest 
power gain. 

Note that the results of Section IIII-BI and Section IIII-CI are 
also valid for the case of MISO RDF. 
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